Percolation thresholds on planar Euclidean relative neighborhood graphs 
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In the presented article, statistical properties regarding the topology and standard percolation on 
relative neighborhood graphs (RNGs) for planar sets of points, considering the Euclidean metric, 
are put under scrutiny. RNGs belong to the family of "proximity graphs", i.e. their edge-set encodes 
proximity information regarding the close neighbors for the terminal nodes of a given edge. Therefore 
they are, e.g., discussed in the context of the construction of backbones for wireless ad- hoc networks 
that guarantee connectedness of all underlying nodes. 

Here, by means of numerical simulations, we determine the asymptotic degree and diameter of 
RNGs and we estimate their bond and site percolation thresholds, which were previously conjectured 
to be nontrivial. We compare the results to regular 2D graphs for which the degree is close to that 
of the RNG. Finally, we deduce the common percolation critical exponents from the RNG data to 
verify that the associated universality class is that of standard 2D percolation. 

PACS numbers: 64.60.ah,64.60.F-,07.05.Tp,64.60.an 
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I. INTRODUCTION 

The pivotal question in standard percolation [J Q is 
that of connectivity. A basic example is 2D random 
bond percolation, where one studies a lattice in which 
a random fraction p of the edges is "occupied". Clus- 
ters composed of adjacent sites joined by occupied edges 
are then analyzed regarding their geometric properties. 
Depending on the fraction p of occupied edges, the geo- 
metric properties of the clusters change, leading from a 
phase with rather small and disconnected clusters to a 
phase, where there is basically one large cluster covering 
the lattice. Therein, the appearance of an infinite, i.e. 
percolating, cluster is described by a second-order phase 
transition. 

There is a wealth of literature on a multitude of vari- 
ants on the above basic percolation problem that model 
all kinds of phenomena, ranging from simple configura- 
tional statistics [27] to "string" -bearing models that also 
involve a high degree of optimization, e.g. describing vor- 
tices in high T c superconductivity 0, |j|, the negative- 
weight percolation problem [3, [3] , and domain wall exci- 
tations in disordered media such as 2D spin glasses 0, H| 
and the 2D solid-on-solid model @. Besides discrete 
lattice models there is also interest in studying contin- 
uum percolation models, where recent studies reported 
on highly precise estimates of critical properties for spa- 
tially extended, randomly oriented and possibly overlap- 
ping objects with various shapes [13]. Further, perco- 
lation phenomena on planar random graphs and their 
duals have been studied extensively in the past [lll - H^ . 
Among the latter graphs are, e.g., 2D Voronoi graphs re- 
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lated to a planar set of points, and their duals, given by 
the Delaunay triangulation of an associated auxiliary set 
of points [lj]- In the latter reference the general inter- 
est of computing percolation thresholds for other random 
systems is declared. 

Here, we consider the Euclidean relative neighborhood 
graph (RNG) for a planar set of, say, N points and de- 
termine the thresholds for bond and site percolation on 
this type of random graph. To prepone some of the de- 
tails given in sect, m note that a graph is subsequently 
referred to as G — (V, E), where V is a set of the nodes 
(represented by a set of N distinct (i-dimensional points, 
called nodes or sites, see Ref. and where E signifies 

the respective edgeset. Considering a particular metric, a 
certain "length" can further be associated with each edge 
(see discussion in sect.Ql]). In a RNG, E contains proxim- 
ity information regarding the close (spatial) neighbors for 
the terminal nodes of a given edge. Here, bond percola- 
tion means that for a given instance of a RNG we occupy 
a fraction p of the graph edges and assess the statistics 
of clusters of adjacent sites connected by occupied edges. 
Examples of bond percolation for an instance of a RNG 
for different values of p are shown in Fig. [Ifa-c). 

The RNG for a given set of points considering an Eu- 
clidean metric was introduced by Toussaint in 1980, see 
Ref. [Hj], who discussed its ability to extract perceptual 
relevant information from a planar set of points. This 
is relevant in the fields of computational geometry and 
pattern recognition, where important questions relate to 
the problem of finding structure behind the pattern dis- 
played by a set of points. RNGs find further applica- 
tion in the construction of planar "virtual backbones" for 
ad-hoc networks (i.e. collections of radio devices without 
fixed underlying infrastructure) , along which information 
can be efficiently transmitted [16H19| . In Toussaints sem- 
inal article it was shown (by means of some illustrative 
examples) that, depending on the precise distribution of 



2 




FIG. 1: Bond percolation on an instance of a relative neigh- 
borhood graph for a planar set of N = 512 points. The sub- 
figures relate to a fraction p of occupied edges, where (a) 
p — 0.72 (sub-critical), (b) p = 0.79 (close to the critical 
point), and (c) p = 0.85 (super-critical). In the subfigures, 
nodes and edges that comprise the largest cluster on the lat- 
tice are colored black, all other nodes and edges are colored 
grey. 

points in the plane, an instance of a RNG might be- 
have similar to a minimum weight spanning tree (MST; 
i.e. a spanning tree in which the sum of Euclidean edge 
weights is minimal, see Ref. 1203) or a Delaunay triangu- 
lation (DT; in a DT, two nodes i, j £ V are joined by 
an edge, if there is a circle passing through them that 
contains no other points k £ V \ {i,j}, see Ref. for 
the underlying set of points. Toussaint showed that for a 
planar set of points, the MST is a (spanning) subgraph 
of the RNG, and further, the RNG is a (spanning) sub- 
graph of the respective DT. Finally, Ref. [15J discusses 
two algorithms that allow to compute the RNG for a 
given set of points, termed ALG-1 and ALG-2. ALG-1 
represents a naive implementation of the RNG construc- 
tion rule (see sect. [II]) that terminates in time 0(N 3 ) 
and is correct for <i-dimensional sets of points as well as 
for non-Euclidean metrics. In contrast to this, ALG-2 is 
rather fast but limited to the planar case and to the Eu- 
clidean metric. Being slightly more "special", ALG-2 is 
based on the observation that in the planar case and for 
an Euclidean metric the RNG is a subgraph of the DT. 
Since there are fast algorithms that allow to compute a 
DT for a planar set of points [ll],[22j (terminating in time 
0(N log(A))), a considerable speedup can be achieved, 
resulting in a worst case running time 0(N 2 ). Further, 
amending ALG-2 by standard techniques to accelerate 
"range-queries" yields an improved worst-case running 
time 0(N log(iV)) (see discussion in sect. HI|) . 

As pointed out above, percolation on the DT of a given 
point set in the plane is well understood. The respec- 
tive thresholds for site and bond percolation can, e.g., be 
found under Ref. [§j|. Regarding the subgraph hierar- 
chy MST C RNG C DT, the question which subgraph 
of the DT still features a nontrivial percolation transi- 
tion was addressed recently, see Ref. [2J] . Intuitively, for 
the MST, the site and bond percolation thresholds are 
1, i.e. the transition points are trivial. Considering the 
RNG for a planar set of points and using the so called 
"method of the rolling ball" , Ref. [24[ established the ex- 
istence of nontrivial site and bond percolation thresholds 
by analytic means. However, numerical estimates for the 



transition points are not provided in the latter reference. 
Here, to elaborate on that, we perform numerical simu- 
lations in order to determine the thresholds of bond and 
site percolation on Euclidean RNGs for planar sets of 
points. Since the RNGs are subgraphs of DTs we can 
expect that the critical exponents that characterize the 
percolation transition on RNGs equal the exponents of 
standard 2D percolation. The transition is hence ex- 
pected to be in the 2D percolation universality class and 
we are primarily interested in the site and bond percola- 
tion thresholds on Euclidean RNGs. 

The remainder of the presented article is organized as 
follows. In section UU we introduce RNGs and the algo- 
rithm we use in order to compute them in more detail. In 
section IIIIl we report on the numerical simulations and 
the analysis of the topological and percolation properties 
of RNGs. Section ITVl concludes with a summary. 

II. MODEL AND ALGORITHM 

RNGs G = (V, E) are based on the concept of rel- 
ative closeness. The nodeset of a A-point RNG is 
given by a set of N distinct d-dimensional points, i.e. 
V = {pi,P2,---,Pn}, where pt = (pi,i, ■ ■ ■ ,Pi,d)- Fur- 
ther, consider a metric L r under which for two points pi 
and pj the distance measure d r (pi,pj) is given by 

d r (pi,Pj) = [V] , \Pi, m ~ Pj.,m\ r ] 1/r ■ (1) 

* — »m= 1 

Then, the edgeset E of the RNG is obtained by connect- 
ing two points pi and pj using an undirected edge {pi, pj} 
if they are at least as close to each other as to any third 
point pfe, see Fig. [DJa). Hence, in order to get joined by 
an edge, the distance d r (pi,pj) of the two points has to 
satisfy the relation 

dij = d r (p t ,pj) < maxldr^p^p^.dripj^pk)} (2) 

for all k — 1 . . . N , k ^ i, j . If Eq. [2] is satisfied, then the 
two nodes are said to be relatively close. In geometrical 
terms, for each pair pi and pj of points, the respective 
distance dij can be used to construct the lune hme(j>i,pj). 
The lune is given by the intersection of two d-dimensional 
hyperspheres with equal radius dij (with respect to the 
prevailing metric), which are centered at pi and pj. If no 
other point pk £ V \ {pi,pj} lies within lxme(pi,pj), i.e. 
if the lune is empty, Eq. [2] holds and pi and pj are thus 
relatively close. In the remainder of the article, if not 
stated explicitly, we consider sets of points in dimension 
d = 2 for the Euclidean metric L^. 

For a planar set of three points (j>\,pi,p^) = (i,j,k), 
the above "selection criterion" for RNG-edges is illus- 
trated in Fig. Ufa). Therein, the individual lunes are 
shown as shaded regions. Since hme(i,j) (encompassed 
by a dashed line) and lune(j, fc) (encompassed by a dot- 
ted line) enclose no further point, the respective pairs 
of nodes are joined by undirected edges. Only lune(i, k) 



3 






(a) 10 



FIG. 2: Illustration of the Euclidean relative neighborhood 
graph (RNG) and its relation to the Denlaunay triangulation 
(DT) and minimum-weight spanning-tree (MST) for the same 
set of points, (a) In a RNG, two points, say i = pi and 
j = pj, are connected by an undirected edge if no third point 
lies within lune(i, j) (see text for more details). Further, a 
slightly less trivial example for N = 30 points, showing (b) 
the RNG (black edges) as spanning subgraph of the DT (black 
and grey edges), and, (c) the MST (black edges) as spanning 
subgraph of the RNG (black and grey edges). 



(encompassed by a solid line) is not empty. It encloses 
the point j, hence the points i and k are not joined 
by an undirected edge. The resulting RNG is thus 
G = (V = {i,j,k},E = So as to fa- 

cilitate intuition, a slightly less trivial example involving 
N = 30 points is shown in Figs.^b-c). In subfigure[2Jb), 
the RNG is highlighted as subgraph of the DT, and, in 
subfigure [2{c) , the MST is indicated as subgraph of the 
RNG. 

Note that a straight-forward implementation of the 
above RNG characteristics can be achieved by consid- 
ering each pair of points (of which there are 0(N 2 )) and 
checking whether one of the remaining N — 2 points lies 
within their lune to rule out that they are RNG neigh- 
bors. This, however, yields an algorithm with running 
time 0(N 3 ), referred to as ALG-1 by Ref. [TpJ (however, 
note that a tremendous speed-up can be achieved by re- 
alizing that already a single point within a given lune 
is sufficient to rule out that the respective lune-defining 
points are RNG-neighbors. Hence, as soon as for the lune 
of a particular pair of points the first such "intruder" is 
identified, one might safely proceed to the next pair of 
points). 

For the planar case and for the Euclidean metric a more 
efficient algorithm, termed ALG-2 (see Ref. [Hj]), can be 
devised. Based on the observation that under the above 
assumptions the RNG is a subgraph of the DT, ALG-2 
can be summarized by the following two steps: (i) con- 
struct the DT Gdt = (Vdt,-Bdt) for the planar set V 
of points, and, (ii) prune the edgeset -Edt of the DT by 
deleting all {pi,pj} e -Edt for which hme(pi,pj) is not 
empty. The latter cleanup phase then results in the edge- 
set E of the RNG for the underlying pointset. So as to 
compute the DT of V in step (i) above, we use the Qhull 
computational geometry library (the DT for a set of 
N points can be computed in time 0(N log(iV)) (2ll,[22T|). 
For the implementation of step (ii) we use the "cell-list" 
method. Therein, the unit square, within which the N 
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FIG. 3: Finite-size scaling analysis for (a) the average running 
time (tpf) of the RNG algorithm for comparatively small sys- 
tem sizes N = 12 . . . 192, (b) the main plot shows the average 
degree (fcjv) of the RNG nodes, and the inset illustrates the 
power-law correction to the scaling behavior of the average 
RNG diameter (R N ). 



points arc distributed uniformly at random, is subdivided 
into Lx L cells (where L = VN), and the 2D cell- indices 
(h,^) = {[L ■ Pi,\\, [L ■ p i:2 \) for all points p { <E V are 
determined. Each cell then is equipped with a list of 
the points it contains. If the lune of a particular pair 
{PiyPj} of points then needs to be checked for empti- 
ness, only a small number mj of cells close to the cells 
with indices (ij.,42) an d (ji, 32) have to be addressed to 
reach all candidate points. Note that the number mj of 
cells to be checked depends on the precise distance dij 
between the respective points. Typically, for large N, 
mj is small for points located in the "bulk" of the unit 
square, and mj can be rather large for points that are 
located along the circumference of the convex hull of V. 
The speed-up achieved by the cell-method is quite im- 
pressive, yielding an improved algorithm ALG-2-CELL 
that terminates in time 0(N log (AT)). This is illustrated 
in Fig.[3][a), where for small systems of size N < 200 the 
average running time (tjv) (averaged over 500 ./V-point 
instances) is shown. The solid line indicated in the fig- 
ure is a fit to a function of the form (t^) oc Nlog(N). 
However, note that for N not too small the data also 
fits well to an effective power-law function that increases 

OCAT1.28W 

Subsequently, we will use ALG-2-CELL to compute the 
Euclidean RNG for 2D sets of points and we compute the 
bond and site percolation thresholds on these graphs. 



III. RESULTS 

In the current section we will use ALG-2-CELL to com- 
pute the Euclidean RNG for planar sets of N = 144(= 
12 2 ) . . . 36864(= 192 2 ) points, where results are averaged 
over 2000 independent RNG instances. In subsect. IHI Al 
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we report on some topological properties of RNGs and 
then we compute the bond and site percolation thresh- 
olds on these graphs. The data analysis for the respec- 
tive bond percolation problem is discussed in detail in 
subsect. IIII B[ and the discussion for the site percolation 
problem in subsect. IIII CI is kept more brief. A visual 
account of bond percolation on an instance of a N = 512 
RNG, i.e. a 2D system of effectively LxL points, is given 
in Fig. QJa-c). 



A. Topological properties of planar RNGs 

First, we discuss two topological characteristics of 
RNGs, namely the average node degree and diameter 
(i.e. the longest among all shortest paths) of the RNG 
in the limit N — > oo. Bear in mind that the nodes of 
a iV-point RNG are distributed uniformly at random in 
the unit square. Thus, one might expect that the scal- 
ing behavior of observables in a 2D RNG depends on the 
effective system length L = N 1 / 2 . 

a. Average degree of the RNG: The scaling behavior 
of the average degree (fcjv) as function of the RNG size N 
is shown in Fig. [3jb) . Therein, the solid line indicates a 
fit to the function (fejv) = fcoo — aN~ b where we find the 
asymptotic average degree koo — 2.5576(3), a = 1.85(5), 
and b = 0.503(5) for a reduced chi-square value Xred = 
0.76 (considering ridof = 6 degrees of freedom). 

b. Average diameter of the RNG: The average di- 
ameter (Rn) of the RNG as function of N is summarized 
in the inset of Fig.[3£b). For a planar graph like the RNG 
one can already expect the approximate scaling behavior 
(Rn) oc N 1 / 2 . Upon analysis we found that the data 
fits best to a function of the form (Rn) = RqN 1 / 2 ^ + 
bN~^^ 2 ], where U) indicates a correction-to-scaling ex- 
ponent. We estimate R = 1.75(2), b = 0.42(1), and 
u = 0.32(3) for a reduced chi-square value x^d = 1-64 
(considering ridof = 6 degrees of freedom). In the in- 
set of Fig. EJb) we aimed to extract the correction to 
scaling according to (Rn)N^ 1 / 2 - R cx N~ u / 2 . The 
numerical value of Rq can also be set into a context 
[Hj]: for the MST on has (R^ ST ) = R^ ST N^ 2 , where 
0.5 < R^ ST < 0.707, see Refs. HHH. Since the RNG is 
a supergraph of the MST one thus might expect to find 



TABLE I: Critical properties that characterize bond and site 
percolation (BP and SP, respectively) on Euclidean RNGs 
for planar sets of points. From left to right: Critical point 
p c (obtained from the analysis of the Binder ratio), critical 
exponents v and f3 obtained from the order parameter, and 
7, obtained from the order parameter fluctuations and the 
scaling behavior of the average size of the finite clusters for 
BP and SP, respectively. 
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RNG-BP 


0.771(2) 


1.33(6) 


0.15(2) 


2.40(6) 


RNG-SP 


0.796(2) 


1.33(6) 


0.14(3) 


2.39(7) 
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FIG 4: Finite-size scaling analysis of the bond percolation 
probability P(p) (i.e. the probability of simultaneous perco- 
lation along both independent directions; see text) on RNGs 
for planar sets of N = 1024 . . . 36864 points, averaged over 
2000 different graph instances. The main plot shows the data 
collapse obtained using Eq.O and the inset illustrates the raw 
data close to the critical point. 
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B. Results for bond percolation on planar RNGs 

To simulate the bond-percolation problem on instances 
of RNGs we implemented the highly efficient, union- 
find based algorithm due to Newman and Ziff, see Ref. 
[27L HH . Therein, initially, each node comprises its own 
(single-node) cluster. We proceed by choosing edges from 
E, one after the other, in random order. For each edge we 
check whether its terminal nodes belong to different clus- 
ters. If this is the case, the respective clusters are merged 
using the "union-by-size" approach. Once all edges are 
considered, the particular RNG instance is completed. 
In principle, this allows to compute observables very ef- 
ficiently with a resolution of 0(1/\E\). For a more clear 
presentation, and so as to be able to compute proper 
errorbars for the observables below, we consider 2000 in- 
dependent RNG instances for a given value of N and 
approximately 80 supporting points on the p axis (in the 
vicinity of the critical point) for which averages are com- 
puted. The observables we consider below can be rescaled 
following a common scaling assumption. Below, this is 
formulated for a general observable y(p, L). This scaling 
assumption states that if the observable obeys scaling, it 
can be rewritten as 

y(p,L) = L- b f[(p- Pc )L^], (3) 

wherein v and b represent dimcnsionless critical expo- 
nents (or ratios thereof, see below), p c signifies the criti- 
cal point, and /[■] denotes an unknown scaling function. 
Following Eq. [31 data curves of the observable y(p, L) 
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recorded at different values of p and L collapse, i.e. fall 
on top of each other, if y(p,L)L b is plotted against the 
combined quantity e = (p — p^)L v and if further the scal- 
ing parameters p Cl v and b that enter Eq. [3] are chosen 
properly. The values of the scaling parameters that yield 
the best data collapse determine the numerical values of 
the critical exponents that govern the scaling behavior 
of the underlying observable y(p,L). In order to obtain 
a data collapse for a given set of data curves we here 
perform a computer assisted scaling analysis, see Refs. 
[23 . [30| . The resulting numerical estimates of the critical 
thresholds and exponents for bond and site percolation 
on planar Euclidean RNGs are listed in Tab. fl] In the 
subsequent paragraphs, we report on the results found 
for different observables. 

c. Percolation probability: In order to provide a 
measure of percolation probability for the 2D RNG, we 
proceed as follows: For each instance of a iV-point RNG 
we first determine the L points that are closest to the 
left, right, top, and bottom boundary. As a sufficient 
condition for percolation along, say, the horizontal di- 
rection, we consider the event that a point on the left 
and the right boundary are part of the same cluster. 
Here, we put under scrutiny the particular event that 
the system simultaneously percolates along both inde- 
pendent directions (other criteria yield similar results). 
The finite-size scaling analysis of the corresponding per- 
colation or "spanning" probability P(p) is summarized 
in Fig. |U Setting b = in Eq. [3] (as appropriate for 
a dimensionless quantity) and restricting the analysis to 
the interval e = [—0.5, 0.5] on the rescaled p-axis, we 
find that p c = 0.771(1) and v = 1.5(2) yield a best data 
collapse with "collapse-quality" S = 0.33 (the numeri- 
cal value of S measures the mean-square distance of the 
data points to the master curve, described by the scaling 
function, in units of the standard error [2!|). Note that 
the numerical value of the correlation length exponent is 
in agreement with the standard 2D percolation exponent 
i/ = 4/3k 1.333. 

d. Order parameter statistics: As a second observ- 
able we consider s max , he. the relative size of the largest 
cluster of points joined by edges. In this regard, a fur- 
ther dimensionless quantity commonly referred to in the 
analysis of phase transitions is the Binder ratio (3lj 
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This ratio of moments scales according to Eq. ((31), where, 
as for the spanning probability above, b — 0. As can be 
seen from the inset of Fig. UKa), the Binder ratio ex- 
hibits a nice common crossing point of the data curves 
for different values of N. The best data collapse (ob- 
tained in the (unsymmetrical) range e £ [—0.1,1.0]) 
yields p c = 0.772(2), and v = 1.4(2) with a quality 
S = 1.22. As evident from the rescaled data (main plot 
of Fig. Eta)), there are rather strong deviations from the 
expected scaling behavior as p < p c . To account for this, 
the scaling analysis is performed in the rather unsym- 
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FIG. 5: Finite-size scaling analyses related to the relative 
S1ZC .S'max of the largest cluster of sites on RNGs for planar 
sets of JV = 1024 . . . 36864 points, averaged over 2000 graph 
instances. The main plots show the data collapse obtained 
according to Eq. [3] and the insets illustrate the raw data 
close to the critical point. The subfigures show different ways 
to analyze s ma x in terms of (a) the Binder ratio b(p), (b) 
the order parameter s ma x(p), and, (c) the fluctuation x(p) = 
iVvar(s max ) of the order parameter. 



metrical interval e £ [—0.1, 1.0] on the rescaled p-axis to 
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accentuate the region p > p c where b(p) seems to scale 
well. As above, the numerical value of the exponent v is 
in agreement with the standard 2D percolation exponent. 
Considering the order parameter 



-Pmax(p) = (Smax(p)), 



(5) 



the best data collapse (obtained in the range e G [—1,1]) 
yields p c = 0.7708(7), v = 1.33(6), and (3 = 0.15(2) with 
a quality S — 0.50, see Fig.[5fb)). If we fix the numerical 
values of the critical exponents to the expected values 
v = 4/3 fa 1.333 and = 5/36 « 0.139 we are left with 
only one adjustable parameter, resulting in the estimate 
p c = 0.7711(6) with S = 0.88. A further critical exponent 
can be estimated from the scaling of the order parameter 
fluctuations xip)i given by 



Xip) = ^[(4ax(p)) " (w(p)> 5 



(6) 



A best data collapse for this observable (attained in 
the range e e [—0.2,0.2]) results in the estimates p c = 
0.769(4), v = 1.30(3), and 7 = 2.40(6) with a quality 
5 = 0.24, see Fig. [SJc). Note that the numerical value 
of the fluctuation exponent is in agreement with the ex- 
pected value 7 = 43/18 2.389. 

e. Average size of the finite clusters: As a last ob- 
servable we consider the average size (S&n(p)) of all finite 
clusters for a particular RNG, averaged over different in- 
stances of RNGs. The respective definition reads pi l32j 



Sfin(p) 



£' s s 2 n s (p) 



(7) 



where n s (p) signifies the probability mass function of 
cluster sizes for a single RNG instance at a given value 
of p. Note that the sums run over finite clusters only 
@jH3 (indicated by the prime), i.e. if the precise con- 
figuration features a system spanning cluster (spanning 
horizontally or vertically or both) , this cluster is excluded 
from the sums that enter Eq. [7] The average size of all fi- 
nite clusters is expected to scale according to Eq.[31 where 
b = —7/^. Restricting the data analysis to the interval 
e G [—1-0, 1.0] on the rescaled p-axis, the optimal scaling 
parameters are found to be p c — 0.770(2), v — 1.36(4), 
and 7 = 2.33(5) with a quality S — 0.61 (not shown). 
Note that here, the numerical values of the extracted ex- 
ponents are in reasonable agreement with the expected 
values and the estimate of the critical threshold for bond 
percolation is consistent with the numerical values found 
above. 



C. Results for site percolation on planar RNGs 

The analysis in terms of the site percolation problem 
was carried out analogous to that of the bond percola- 
tion problem in the preceding subsection. However, we 
here list only the estimates of the critical points obtained 
from the data collapse analysis for the different observ- 
ables. In this regard we have p c = 0.799(1) (percolation 



probability), p c = 0.796(2) (Binder ratio), p c = 0.795(1) 
(order parameter), p c = 0.794(4) (order parameter fluc- 
tuation; see Tab. HJ, p c — 0.798(3) (average size of the fi- 
nite clusters). Further, the critical exponents v = 1.33(6) 
and P — 0.14(3) obtained from the order parameter and 
7 = 2.39(7), obtained from the scaling behavior of the 
average size of the finite clusters, are listed in Tab. [I] 



IV. DISCUSSION AND SUMMARY 

In the presented article we have closely investigated 
the statistical and percolation properties of planar Eu- 
clidean RNGs via numerical simulations. In regard of the 
subgraph hierarchy MST C RNG C DT, recently the 
question was raised which subgraph of the DT (on which 
the percolation problem is well studied) still features a 
nontrivial percolation transition [24| . Intuitively, for the 
MST this is not the case. For the RNG previous analytic 
studies established the existence of nontrivial site and 
bond percolation thresholds 24], but no numerical esti- 
mates where provided. Here we quote p^bond = 0.771(2) 
and p™^ — 0.796(2), obtained by means of finite-size 
scaling analyses in terms of the "data-collapse" tech- 
nique. Further, we also deduced the critical exponents 
that govern both percolation transitions on the RNG and 
found them to be consistent with those that describe the 
standard 2D percolation phenomenon (as expected). So 
as to yield maximally justifiable results through numeri- 
cal redundancy, we considered various observables to es- 
timate the critical points and exponents. 

As discussed in subsect. MI A| the asymptotic aver- 
age degree for the RNG reads fcoc = 2.5576(3). In or- 
der to put the above critical points into a context, we 
might attempt to compare them to the threshold val- 
ues for regular 2D lattices with a similar degree. E.g., 
the site and bond percolation thresholds for the (3,12 2 )- 
Archimedean lattice, having degree k = 3, read p c ,bond = 
0.74042081(10) H and p c , sito = 0.807900764.'.. [34|. 
Both are actually not thus far from the respective thresh- 
olds on the RNG. (Further, for the k = 3 Martini lat- 
tice one has Pn,bond = 0.707107 . . . [35J and p CjS itc = 
0.764826... |3g|.) Regarding 2D random lattices, site 
and bond percolation on the Voronoi tesselation of a 
planar pointset, also having degree k = 3, give rise to 
the thresho ld v alues p c ,bond = 0.666931(5) and 
0.71410(2) 12]. In addition, site percolation on planar 
$ 3 random graphs result in p c , S it c = 0.7360(5) [13 ]. In 
comparison, the estimates from the latter random graphs 
are less close to the estimates for the RNG as compared 
to the (3,12 2 )-Archimedean lattice thresholds. 

As pointed out in the introduction, RNGs are discussed 
in the context of the construction of planar "virtual back- 
bones" for ad-hoc networks that guarantee connectedness 
of all considered nodes [l6l - [l9j . In this regard, from a 
point of view of stability, it would be interesting to quan- 
tify how susceptible RNGs are with respect to a random 
failure and targeted removal of nodes and to compare 
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the results to other "proximity graph" instances which 
are discussed in the same context. Such investigations 
are currently underway 37]. 
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